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ABSTRACT 

We compute the structure and degree of neutronization of general relativistic mag- 
netohydrodynamic (GRMHD) outflows originating from the inner region of neutrino- 
cooled disks. We consider both, outflows expelled from a hydrostatic disk corona and 
outflows driven by disk turbulence. We show that in outflows driven thermally from a 
static disk the electron fraction quickly evolves to its equilibrium value which is domi- 
nated by neutrino capture. Those outflows are generally proton rich and, under certain 
conditions, can be magnetically dominated. They may also provide sites for effective 
production of ^^Ni. Centrifugally driven outflows and outflows driven by disk turbu- 
lence, on the other hand, can preserve the large in-disk neutron excess. Those outflows 
are, quite generally, subrelativistic by virtue of the large mass flux driven by the addi- 
tional forces. 

Subject headings: accretion, accretion disks - gamma rays: bursts - MHD - nuclear 
reactions, nucleosynthesis, abundances - relativity 



1. Introduction 

The structure and nuclear composition of outflows from a disk accreting rapidly onto a central 
black hole is a topic of considerable interest. Such systems are believed to form following the 
collapse of rotating massive stars in failed supernovae events (coUapsars), or the merger of neutron 
stars and black holes, and there is some evidence that link them to GRBs. For the range of accretion 
rates {Mace ~ 0.01 — 10 Mq yr~^) and viscosity parameters (c^vis ~ 0.01 — 0.1) anticipated, the disk 
surrounding the black hole is sustained at MeV temperatures and its midplane density is high, in 
excess of lO^'' gr cm~^ (e.g., Popham et al., 1999). Under such conditions the inner parts of the 
disk, within about 10 to lOOr^ depending on parameters, consist of neutron rich matter that cools 
predominantly via emission of MeV neutrinos (Popham et al. 1999; Di Matteo et al. 2002; Kohri 
& Mineshing 2002; Pruet at al. 2003; Chen & Beloborodov 2007, hereafter CB07). For M^cc ~ 0.1 
and Ovis ~ 0.03 the neutron-to-proton ratio within lOr^ exceeds 30 (CB07). 

Viscous and neutrino heating of the upper disk layers drives a powerful wind from the disk. 
This wind is most likely baryon rich and expands at sub or mildly relativistic speeds except, perhaps. 
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inside a core containing the putative baryon poor GRB jet. The baryon rich wind emanating from 
the disk is a potential site for efficient nucleosynthesis (Pruet et al. 2003, 2004), and may account 
for the SN-like features observed in several long GRBs. It may also play an important role in the 
collimation of the central GRB-producing jet (Levinson &: Eichler 2000; Alloy et al. 2000; Rosswog 
&; Ramirez- Ruiz 2003; Alloy et al. 2005; Bromberg & Levinson 2007). The latter may be launched 
magnetically from the very inner regions of the disk (e.g., Vlahakis & Konigl 2003; Levinson 2006) 
or may be produced via a Blandford-Znajek mechanism (Levinson &; Eichler 1993; Meszaros & Rees 
1997; Van Putten 2001; Van Putten & Levinson 2003). If, as often assumed, the GRB jet picks 
up nuclei from the disk (e.g., Derishev et al. 1999; Beloborodov 2003) then it may contain matter 
with a large neutron excess that is likely to affect its dynamics (Fuller et al. 2000; Vlahakis, Peng 
&; Konigl 2003; Rossi et al. 2005), as well as some of the characteristics of the afterglow emission 
(e.g., Derishev et al. 1999; Bachah & Meszaros 2000). 

Whether disk outflows can preserve the large in-disk neutron excess is yet an open issue. 
Pruet et al. (2004) constructed a simple model of hydrodynamic disk winds and argued that in 
such winds the electron fraction typically evolves to 1^ ^ 0.5. They also pointed out that lower 
values of the asymptotic electron fraction are expected in centrifugally driven winds. Their model 
assumes effectively that the outflow is ejected from a hydrostatic corona, where the density is high 
and the flow velocity is very small. As proposed by Beloborodov (2003), turbulent mixing can lift 
up neutron-rich matter from the inner disk layers to its surface on a timescale much shorter than 
the neutronization time. A fraction of this circulating material may then be ejected as an outflow, 
in which case the conditions at the base of the flow may be vastly different than those expected in 
the case of an outflow from a hydrostatic corona. 

In a recent paper (Levinson 2006; hereafter LE06) a model for GRMHD outflow from a 
neutrino-cooled disk has been developed and employed to calculate the structure of the outflow 
in the sub-slow magnetosonic region and the mass loading of the outflow. The method used to 
compute the mass flux is outlined in f|2] below. The principle conclusion drawn in that paper is that 
ejection of relativistic outflows from the innermost disk radii, within several or so, is possible 
in principle for certain magnetic field configurations even in non-rotating black holes, provided 
the neutrino luminosity is sufficiently low, 10^^ erg s^^, and the magnetic field is sufficiently 

strong, B ~ 10^^ G. The conditions found to be optimal for the launching of an ultra-relativistic jet 
are also the conditions favorable for large neutron-to-proton ratio in the disk. However, the compo- 
sition profile of the outflow was not computed in LE06, and the question whether the outflow can 
retain a large degree of neutronization was left open. In this paper we incorporate the evolution 
equation for the electron fraction into the model developed in LE06, and compute numerically the 
structure and nuclear composition of the outflow beneath the slow magnetosonic point for a range of 
conditions in the disk. We specifically consider outflows emanating from a hydrostatic disk corona 
and steady outflows driven by disk turbulence. 
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2. 



GRMHD Disk Outflow Model 



We consider a stationary, axisymmetric MHD wind expelled from the surface of a hot, mag- 
netized disk accreting onto a Schwarzschild black hole. The range of conditions in the disk is 
envisaged to be similar to that computed by CB07 for a black hole of mass SMq, accretion rates 
Mace ~ 10^^ — 1 Mq s^^ and viscosity parameters Qvis = 0.01 — 0.1. Those authors generalized ear- 
lier work by Popham et. al. (1999) and Pruet et al. (2003) to incorporate calculations of electron 
degeneracy and nuclear composition self-consistently, which affect the conditions in and neutrino 
emission from the disk. Under the conditions envisioned above the dominant cooling mechanism 
in the disk is neutrino emission, with a total luminosity of Ly + Ly ~ 0.04MaccC^ for accretion 
rates above the ignition values (CB07). The neutrino luminosity is considerably higher if the black 
hole is rapidly rotating. The major fraction of the neutrino luminosity is generated in the inner 
disk regions, within 10 Schwarzschild radii or so, although for very low ctviss the neutrino emission 
may extend to much larger radii. The matter in the inner disk regions is typically neutron rich. 
The electron fraction decreases, quite generally, with increasing Mace and decreasing Ovis- For 
example, for Mace = 0.2 Mq s~^ the electron fraction changes from Yf, ~ 0.03 to Yf. ~ 0.15 as QvIs 



The model outlined in LE06 calculates the structure of the GRMHD outflow below the slow 
magnetosonic point for a given magnetic field geometry, treating the neutrinos emitted from the 
disk as an external energy and momentum source. To simplify the analysis the neutrino source has 
been taken to be spherical with radius R^. The model is characterized by three parameters: the 
black hole mass Mbh, the neutrino flux L^/R"^, and the mean neutrino energy < Ey >. It solves 
a set of coupled ODEs that are derived from the general relativistic energy- momentum equations, 
describing the change along a given streamline, ^{r,0) = const, of the specific energy E , specific 
entropy s, and poloidal flow velocity Up-. 



Here ' denotes the derivative along streamlines, u°'da, denotes the source terms associated with 
energy and momentum exchange with the external neutrino source, u" is the outflow 4-velocity, p is 
the baryon rest mas density and T is the temperature. The denominator on the right hand side of 
eq. ([3]) is given explicitly as D = — (a^ — — M^)^(Up— ■u|j^,^)(Mp — M|.^^)/tt^, where ua, usm and 
upM are the Alfven, slow and fast magnetosonic wave speeds, respectively, a is the lapse function, Q 
is the angular velocity defined below, R is the cylindrical radius and M = Up/uA is the Alfven Mach 
number. The term F can be expressed as F = Ci(lni?p)' + C2(lna)' + C3(lni?)' + C4(ln£')' + C5(lns)', 



is varied from 0.01 and 0.1 (CB07). 



2.1. Basic Flow Equations 



{p/mN)kBTs' = -Uaq", 
pc^E' = -qt, 

(In Up) = — . 



(1) 
(2) 

(3) 
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where the coefficients Cfc are functions of the flow parameters, given exphcitly in LE06, and Bp is 
the poloidal field component. Since the derivatives {In Bp)' and (Ini?)' depend on the magnetic 
field geometry which is unknown a priori, additional equation is needed. Our approach is to invoke 
a given field geometry. To examine the dependence of mass flux on the latter, we obtained solutions 
for different magnetic field configurations, focusing particularly on split monopole and r self-similar 
geometries. Equations ([I])-(l3]) are augmented by an equation of state for the mixed fluid of baryons, 
photons and electron-positron pairs. In addition there are three integrals of motion of the MHD 
system: the mass-to-magnetic flux ratio ?/(^'), the angular velocity of magnetic field lines $^(^'), 
and the specific angular momentum L{^). The two invariants ^{^) and L(^) are fixed by a choice 
of boundary conditions. The mass flux r](^) is an eigenvalue of the problem, and is determined by 
the regularity condition at the slow magnetosonic point. 



2.2. The Electron Fraction 

For the range of conditions considered below the matter in the sub-slow magnetosonic region 
consists mainly of free nucleons. The total baryon density is then given hy p = rUpUp + mnUn, 
where n„ and Up denote the number density of neutrons and protons, respectively, and m„, rup the 
corresponding masses. The neutron-to-proton ratio is related to the electron fraction Yg through 

W/^) = ^^, (4) 

Tip -\- Tlji 

and is determined by a competition between the following reactions: 

e~+p^n + Ve, (5) 
e^ + n^p + Ue- (6) 

Neutron decay is negligible due to the long life-time of the neutron compared with the characteristic 
timescales involves. Lepton capture on heavy nuclei can be ignored since as stated above (and will 
be confirmed below) in the regime considered here the matter consists mainly of free nucleons. 
The reactions Ve + T^e e"*" + e~ are typically unimportant at the characteristic densities and 
temperatures involved, nonetheless, they are incorporated for completeness in the source terms 
that appear on the right-hand side of eqs. ([ID and d^D (see LE06 for details). 

The change of the electron fraction along a given streamline is determined by 

where Ag-p, Ae+„, App, X^n-, are the rates for the forward and reverse reactions in eqs. ([5]) and ([6]), 
and ' denotes again the derivative m"9q,. 

The rates for electron and positron capture are given by (Fuller et al. 1980), 

/>oo 

\^-p = k w'^{w-A/mefG^{l,w)S^{l-Su)dw, (8) 

J A /me 
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/oo 
w'^{w + /^/mefG+{l,w)S+{l- Sv)dw (9) 

where k ~ 6.414 X 10~^ is comparative half-life related to the Gamow- Teller and Fermi matrix 
elements, A = m„ — m,p = 1.293 McV is the neutron-proton mass difference, vrie is the electron 
mass, w is the total energy (rest mass and kinetic energy) in units of 'm,,,c^, and 5"+ arc the 
electron and positron distribution functions, respectively, and G±{Z = l,w) are the Coulomb 
correction factors, given in Fuller et al. (1980). For the range of temperatures and densities 
considered the e^ pairs and photons are in perfect thermodynamic equilibrium in the sub-slow 
magnetosonic region, owing to the huge optical depth there. The functions S± are then given by 
S±{T, jl) = [expK^y lb lJ-)/0} + 1]^"*^, where d = /m(,(? and ^ is the electron chemical potential 
(measured in units of mgC^). The latter is determined from the charge neutrality condition: 

-^Ye = n_ - n+ = ^ (^)' r [S_(T,/x) - S+(r, /x)]/dp, (10) 

with p = [viP' — 1)^/^ being the electron/positron momentum in units of rWeC. 

As mentioned above, we assume for simplicity that the neutrino source is spherical with some 
characteristic luminosity L^,, mean neutrino energy < E,^ >, and radius Ri,. In cases where the disk 
is z/-opaque Ri, denotes the radius of the corresponding z^-sphere. Due to the potentially different 
opacities for and absorption the luminosity Lp, mean energy < Ep > and photosphcric radius 
Rp of the antineutrinos may differ from that of neutrinos. The details depend on the temperature 
profile in the neutrino production zone, which is governed by the dissipation mechanism in the 
disk. In z^-transparent disks we naively anticipate Lj^/Lp =< E,j > / < Ep >= R^/Rp = 1. The 
rates for the reactions + n ^ p + e~ and i/g + p — n + e+ at the base of the outflow can be 
approximated as (Qian and Woosley 1996) 

A„ = i±j^G||(.„ + 2A + ^)(l-..), (11) 

where a ~ 1.26, = 5.29 x 10~^^ cm^ MeV~^ is the Fermi coupling constant, x^(^p) = (1 — 
i^^^P^/r^)-^/^ is a geometrical factor, and ej,(^) =< E"^^-^ > / < E^^^p^ >. Detailed calculations of 
the neutrino spectrum emitted by the disk is beyond the scope of this paper. In what follows we 
shall assume for simplicity that the z^e and i?e are emitted from the same region with the same 
luminosity and spectrum, that is, we take Ri, = Rp, Ly = Lp and < Ey >=< Ep >. The energy 
moments < Ey > and depend on the shape of the neutrino spectrum. For a blackbody spectrum 
< Ey >= 3.15kBTy, = AAksTy, where Ty denotes temperature of the neutrino source, whereas 
for a neutrino transparent source < Ei, >~ bksTiy, ~ GkBT^ (Beloborodov, 2003). The free 
parameters L^, T^, and R^ are estimated by employing the results of CB07. 
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2.3. Conditions at the Flow Injection Point 

In LE06 it has been assumed that the outflow connects to a hydrostatic disk corona where the 
density is high and the entropy per baryon, 

s = 8.7 + HTlily/p^) + 0.53(r|,,v'/P9), (13) 

is relatively small. Under this assumption, the integration of eqs. ([I])-® starts at sufficiently dense 
disk layers where the light fluid pressure, pi = (ll/12)ar^, roughly equals the baryonic pressure, 
Ph = pkBT/rUp. The mass flux 77 is then adjusted iteratively by changing the boundary value of the 
poloidal velocity UpQ, until a smooth transition across the slow magnetosonic point is obtained. 

Different boundary conditions may apply if the disk is turbulent. Beloborodov (2003) and 
Pruet et al. (2003) proposed that turbulent mixing may quickly lift up neutron-rich matter from 
the inner disk layers to its surface, and that a small portion of circulating material will then be 
picked up by the outflow. Since vertical mixing is expected to occur over the sound crossing time, 
tjnix — H/cs ~ where H is the disk scale height, which for disk temperatures < 8 MeV 
is smaller than both, the neutronization time and neutrino cooling time (Beloborodov 2003), we 
anticipate the temperature and electron fraction at the outflow injection point to be roughly equal 
to their values at the disk midplane. The disk turbulence is expected to be subsonic and, therefore, 
the poloidal velocity Upo of circulating matter which is ejected into the flow is anticipated to be a 
fraction of the sound speed inside the disk, viz., Upo ^ Cs ~ O.OScTmcV- This is typically below the 
slow magnetosonic speed at the flow critical point (as will be confirmed below). Consequently, under 
the assumption that the outflow above the turbulent layer is steady (averaged over times longer 
than the mixing time) it should pass through the slow magnetosonic point. This again determines 
the mass flux, as in the case of outflows that connect to a hydrostatic corona. However, in this 
scenario the outflow starts from a layer of much smaller density and much higher specific entropy 
than in the solutions obtained in LE06. To account for such cases, we relax the constraint imposed 
in LE06 on the specific entropy at the fiow injection point. Specifically, for any given choice of 
initial temperature, Tq = T{r = tq), and electron fraction, Y^o = l^(r = tq), we construct a family 
of transonic solutions that are characterized by one parameter: the initial entropy sq = s{Tq,pq) 
(or equivalently the initial density pq). 



3. Analytic Treatment 

In regions where the outflow is radiation dominated, viz., pi/pb > 1, the density satisfles 

P9 < 0.13r|,,^. (14) 

The degeneracy condition pg > 0.2T^j^y implies that in this region of the flow the electrons are 
non-degenerate. Under such conditions the capture rates for electrons and positrons are given 
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approximately by X^-p — Ag+n — 0-45T^gy s ^ provided TmcV ~ 1 (Qian & Woosley 1996). Using 
eqs. (jlip and (jl2p we then obtain 

^ ^ 2 X 10^ 6.,Mey R-J (1 - x.). (15) 



Ae~p 



The latter ratio can also be expressed in terms of the neutrino heating rate, q = 5Lj^52e^ MeV-^fd^ ~ 
Xy) MeV s^"*^ baryon"-*^, as 



Ag+„ \ 1 MeV s ^ baryon 

For accretion rates above the ignition value Li, = {Ly + Ly)j2 ~ 0.02MaccC^ in case of accretion 
onto a 3M0 Schwarzchild black hole (CB07). Adopting for illustration Mace = 0.2Mq s~^ and 
Rv% = 5 yields Li,^2/R^q ~ 0.03. For the latter choice of Mace CB07 obtained a peak temperature 
of Ty = 2.5 MeV for Oviss ~ 0.03 and Ty = 2> MeV for Oviss = 0.1. With e^^MeV = ^Ty^MeV we then 
have at the base of the flow A^n/Ae+„ = 0.75 and 0.35, for Oviss ~ 0.03 and Oviss = 0.1, respectively. 
However, as the flow decelerates the temperature drops and this ratio quickly rises. For the range 
of disk parameters considered here, Ly^2 ~ 0.2, e^^MeV > 10, Tq^meV ~ 3, and R^q ~ 3, we find, 
quite generally, that the asymptotic electron fraction is governed by neutrino capture. 

Denoting by dl the proper length along a streamline in units of Ry, and neglecting capture, 
we can write eq. ([7]) as 

dY, ^ KnRy^^ _ _ W^Ye. (17) 

dl Up Up 
Adopting ey = bksTy and using eq. pT]) we obtain at the flow injection point, 

~ 0.02c L,52 Ty^MeV Ryi u-^ , (18) 

Up 

where c is the speed of light. From eq. (jl7p it is evident that the asymptotic value of the electron 
fraction will remain near its initial value provided the outflow time, ~ Ry/upo, is shorter than 
the neutronization time, A~^, which requires 

Upo > 0mcLy52Ty^MevRyi- (19) 

In the following, it is shown that this condition is satisfied in the case of centrifugally driven winds 
and outflows driven by disk turbulence. In thermally driven winds the outflow time is much longer 
than that implied by condition (|19p . The electron fraction then quickly evolves to its equilibrium 
value, Ye^eq — Xyn/{Xvn + App). For a symmetric z^e and t'e emission eqs. ([TT]) and (fT2|) yield 
Xyn > Xyp due to the threshold effect, resulting in asymptotic proton excess (l^,oo > 0.5) in those 
outflows. Our detailed calculations yield Ve.oo ^ 0.5 also in cases where capture dominates the 
evolution of the electron fraction in the wind. 

At the slow magnetosonic point the sound speed age is equal to a modified escape speed. 
Typically a^c ~ 0.1c (LE06). Since ~ Api/Spc^ in the subslow region it implies that the density 
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at the injection point must satisfy po ~ lO^^Mey §^ ^- This, in turn, imphes a (one sided) mass 
flux of 

M ~ vrr^pupo ^ lo30r|^^^(npo/0.1 c) gr cm'^ (20) 

for an injection radius tq — Ry = 3 x 10^ cm. Those rough estimates are confirmed by numerical 
integration of the full set of equations. The above considerations suggest that for explosion energies 
inferred in GRBs neutron rich outflows can only have modest Lorentz factors (F ~ a few), unless 
somehow the turbulent mixing process manage to selectively accelerate only a very small fraction 
of the neutron rich matter that is lifted to the disk surface to velocities well in excess of the escape 
velocity. 

4. Numerical Results 

Equations ([I])-©, ([7|)-(fT2|) have been integrated numerically, as in LE06, using a split monopole 
magnetic field geometry. Each field line in this configuration is characterized by two parameters: 
the inclination angle 9 of the field line, here measured with respect to the symmetry axis, and the 
radius Rq at which the field line intersects the equatorial plane. In the following we consider, in 
turn, outflows from a steady disk and outflows driven by disk turbulence. 

4.1. Steady Outflow from a Hydrostatic Disk Corona 

The integration of these models starts in the dense disk layers, where the pressure is dominated 
by the baryons. The initial temperature Tq is taken to be equal to the temperature of the neutrino 
source, viz., Tq = Ty. For a given choice of Ty the density at the origin, /jq, is chosen such that 
the light fluid pressure, pi = (ll/12)aTQ , does not exceed the baryonic pressure, = poksTQ/mp. 
The initial entropy sq = s(To,po) is then determined using eq. p^ . To verify that the solution is 
insensitive to our choice of boundary conditions, each integration has been repeated several times 
for a given choice of our model parameters, each time with a different value of po- We find that as 
long as pi{Tq) < p;,(To,po) the results are indeed highly insensitive to our choice of po, except for 
the initial value of the electron fraction which at these densities equals its equilibrium value. For 
our fiducial model we choose the electron fraction at the injection point to be l^o = 0.1, which is 
the equilibrium value at a density p ~ lO^'^ gr cm~^ and temperature T ~ 2 MeV (Beloborodov, 
2003). 

In general there are two distinct regimes. The regime of unstable equilibrium corresponds 
to magnetic field lines having inclination angles 9 > 7r/6. Along such field lines the outflow is 
centrifugally driven and can be initiated even in the cold fluid limit (Blandford &: Payne 1982; see 
LE06 for a generalization of this result to the general relativistic case). The slow magnetosonic 
point in this case is located very close to the disk surface, at « ^o^ where the density is high. 
Because of the high critical density the mass flux along such field lines is large and depends only 
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weakly on the neutrino luminosity emitted from the disk. For reasonable magnetic field strengths 
those outflows are typically subrelativistic. The regime of stable equilibrium corresponds to field 
line inclination angles 6 < vr/G. Along such field lines the mass flux is thermally driven, similar 
to the case of a spherical wind, and is a sensitive function of the neutrino luminosity. The slow 
magnetosonic point in this regime is located higher above the disk, at heights Zsm ^ (see fig. 
1). The heating of the wind by the escaping neutrinos results in a steep rise of the entropy per 
baryon, s, during the initial acceleration phase, after which it saturates. The asymptotic value of s 
is larger for lower L^/R^, but does not seem to reach extreme values {s ^ 100). Adiabatic cooling 
is important along field lines of sufficiently small 6, rendering the critical density and temperature 
substantially smaller than in the case of centrifugally driven winds. As a result the mass flux is 
strongly suppressed. For 9 < 15° the outflow is magnetically dominated if the strength of the 
poloidal magnetic fleld at the disk surface is Bp ~ 10^^ Gauss, as often envisioned. The magnetic 
energy per baryon depends sensitively on the inclination angle of the fleld line near the surface and 
can reach very large values (LE06). 

Examples of thermally driven winds are exhibited in flg. 1, where the flow parameters are 
plotted as functions of the height above the disk midplane, for tan^ = 0.2 and vq = = 3rg. The 
left panel corresponds to a disk temperature Tq = 2 MeV and neutrino luminosity = 6x 10^^ erg 
s~^, and the right panel to a disk temperature Tq = 3 MeV and neutrino luminosity = 2 x 10^^ 
erg s~^. The initial density and electron fraction in both examples are po = 10^^ gr cm^'^ and 
YeQ = 0.1, respectively. The disk is transparent in both cases, and so the average neutrino energy 
is taken to be < Ei^ >~ S/cTq. The mass flux, as defined in eq. ()20p . is found to be about 10^^ gr 

in the former case and 10^'^ gr s~^ in the latter, corresponding to a net energy per baryon of 
E = 500 and E = 5, respectively, for a poloidal magnetic field strength of BpQ = 10^^ Gauss. As 
stated above, the solution is highly insensitive to the choice of initial density po, provided pi/pb < 1- 
Changing the latter from 10^^ to 10^ gr cm~^ altered the mass flux and asymptotic 1^ by less than 
one percent. It also had little effect on the profiles of the flow quantities. 

From flg. 1 it is seen that the electron fraction reaches its final equilibrium value close to 
the base of the flow. At the origin, the rates for neutrino capture and e^ capture are comparable 
and Ye adjusts instantaneously to its local equilibrium value, l^^eq(/007 ^o) — 0.3. As the density 
and temperature drop the ratio Xun/^e-p increases, leading to a further increase of Ye until it 
reaches its flnal equilibrium value ^6,00 — ^un/i^un + Ap) ^ 0.5 . This rapid evolution of Yg is 
characteristic to all thermally driven outflows, and is a consequence of the fact that the flow time, 
tf = {tq/c) J dl/up ~ RyjupQ^ implied by the critical condition is considerably longer than the 
neutronization timescale A^^. In the examples depicted in fig. 1 the poloidal velocity at the fiow 
injection point is UpQ = 4.5 x 10^^ (left panel) and UpQ = 4.5 x 10"^ (right panel), well below the 
value required to retain low Ye (see eq. jl9j). The asymptotic electron fraction depends to some 
extent on the disk temperature T^, owing to the threshold effect (i.e., the dependence of the rates 
Xun and Xpp in eqs. [TT] and [12] on the neutron-proton mass difference A). In general, l^,eq tends 
to approach 0.5 as T^, increases. It should be kept in mind, however, that the uncertainty in the 
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luminosity and mean energy ratios of v and v should be reflected in the final value of Ye- 

As stated in ^2.21 above, our analysis ignores the presence of heavy nuclei. As a consistency 
check, we plotted in fig. 1 the equilibrium value of the free nucleon mass fraction, X^-^^^i^^ p), taken 
from Woosley & Baron (1992) (Xfj-ee = 1 implies free nucleons only). This estimate of Xfj-ee is 
accurate at the base of the flow, where the matter is in rough kinetic equilibrium. Further up, 
where the flow velocity exceeds values at which the expansion time becomes shorter than the weak 
interaction timescale it underestimates the fraction of free nucleons. As seen, our neglect of lepton 
capture on heavy nuclei in the relevant phase of the flow, where the change in the electron fraction 
is signiflcant, is justifled. 

Example of a centrifugally driven outflow is shown in flg. 2, for two different inclination angles. 
As seen, the electron fraction evolves rather slowly. This is a consequence of the short expansion 
time in the sub-slow region. In general, the asymptotic electron fraction would depend on the 
structure of the flow above the slow magnetosonic point. However, given the high poloidal velocity 
at the critical point {upc ~ 0.04c in both cases shown in flg. 1) we expect only little evolution above 
the slow point. The general conclusion is that centrifugally driven outflows can retain low values 
of Ye- This trend is consistent with that found by Pruet et al. (2004). The angular velocity of 
magnetic field lines in those calculations is O = 0.95r2ii', where = {'"'g/'^ToY^'^ is the Keplerian 
angular velocity at the disk midplane. The corresponding mass fiuxes are M ~ 0.1 and 0.2 Mq s~^, 
for 9 = 30° and 9 = 60°, respectively. Such high mass loss rates imply that a significant fraction 
of the accreted matter will be expelled from the disk before reaching the black hole, which should 
affect the disk structure considerably. In particular, the disk is likely to become sub-Keplerian, and 
this would feed back on the flow. As illustrated in LE06, the mass flux is highly sensitive to the 
choice of 0. For example, repeating the calculations in fig.l with i7 = 0.85 yielded a mass flux 
smaller by about an order of magnitude. Interestingly, the location of the slow magnetosonic point 
was not altered at all, and the poloidal velocity at the critical point changed by only 50%. The 
reason is that the outfiow accelerates faster. As a consequence the asymptotic electron fraction has 
not changed significantly. 

4.2. Steady Outflow from a Turbulent Disk 

As explained above, this model assumes that turbulent mixing lifts neutron rich matter from 
the inner disk layers to the surface on a timescale much shorter than the neutronization time. The 
temperature Tq and electron fraction YeQ at the flow injection point are then taken to be equal 
to their values at the disk midplane, as before. However, the entropy per baryon at the injection 
point is allowed in this model to be much higher than in the case of outflows that emerge from a 
hydrostatic corona, and is treated essentially as a free parameter. The only restriction imposed on 
the solution is that it should start sub-slow, that is, the slow magnetosonic Mach number at the 
flow injection point must be be smaller than unity. 
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In a turbulent disk the magnetic field is likely to be disordered. In that case the stream function 
should be considered as defining the direction of streamlines. We stress that our calculations are 
applicable also to hydrodynamic flows. It is worth noting that a flow can, in principle, be Poynting 
flux dominated even in case of disordered magnetic fields (e.g., Proga et al. 2003), although, 
as discussed below, for the outflows considered here the magnetic energy per baryon in typically 
smaller than unity. 

Sample results are shown in fig. 3. As expected, the evolution of the electron fraction is 
suppressed as the initial Mach number of the injected flow is increased. We find neutron excess 
(y^oo < 0.5) in solutions for which the slow magnetosonic Mach number at the flow injection point, 
MgM, exceeds 0.1 roughly. At the same time the mass flux also increases. We find that solutions 
that retain Yg < 0.5 have a relatively large mass flux, M ^ 10^*^ gr s~^, in accord with eq. ()20p . We 
therefore conclude that neutron rich outflows expelled from a turbulent disk are likely to be sub or 
at best mildly relativistic. Note that neutrino heating is unimportant (negligible change in specific 
entropy in all the solutions exhibited in fig. 3). This is because of the relatively short expansion 
time at the base of the fiow, tfX^n < 1- Thus, the structure of those outflows depends on the disk 
temperature, but is independent essentially of the neutrino luminosity ( does depend on L^, of 
course) . 



5. Summary and Conclusions 

We have computed the structure and the neutron-to-proton ratio in GRMHD outflows from 
a neutrino-cooled disk accreting onto a Schwarzschild black hole. We considered both, outflows 
emanating from a hydrostatic disk corona and steady outflows driven by disk turbulence. The 
main results are: 

1. In thermally driven outflows that emerge from a steady disk the neutronization time at the 
base of the outflow is much shorter then the outflow time. In this case the electron fraction quickly 
evolves to its equilibrium value which, for the parameter regime explored above, is dominated 
by neutrino capture. For z^-transparent disks this implies asymptotic proton excess {Ye > 0.5). 
The mass flux driven by neutrino and viscous heating along magnetic fleld lines inclined at small 
angles to the vertical {6 ^ 15°) is found to be rather small for low viscosity disks (oviss ~ 0.03) and 
moderate accretion rates (Mace ~ 0.1). If the outflow energy is extracted magnetically from the 
disk with a luminosity comparable to that observed in GRBs, then it can in principle accelerate to 
a high Lorentz factor. Luminous disks (L^, > 10^^ erg s~^, ,^ 3 MeV) give rise to a much larger 
mass loss rate (M ^ 10~^ Mq s~^). Such outflows can in principle synthesize the amounts of ^^Ni 
required to explain the lightcurves of the associated SN (Pruet et al. 2004). 

2. Centrifugally driven winds can retain large neutron excess by virtue of the much shorter 
expansion time. In those outflows the slow magnetosonic point occurs very close to the disk surface, 
where the density is high. The large critical density results in a substantial baryon loading of the 
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flow, rendering it subrelativistic. The large mass loss rates obtained imply that those outflows 
should considerably affect the disk structure (this conclusion may hold true in general if magnetic 
extraction of angular momentum is significant). A complete treatment requires self-consistent 
solutions of the disk and the outflow. 

3. Outflows driven by disk turbulence can also retain large neutron excess if the slow magne- 
tosonic Mach number at the origin satisfies Msm ~ 0.1. In this scenario a fraction the circulating 
material is ejected into the flow over a timescale much shorter than the neutronization time (Be- 
loborodov 2003; Pruet et al. 2003). The in-disk temperature and composition are then expected to 
be preserved at the base of the flow. As a result, the asymptotic entropy per baryon is larger than 
in outflows from a steady disk (see fig. 3). The poloidal velocity at the flow injection point is likely 
to be a fraction of the sound speed inside the disk, which is typically large enough to suppress the 
evolution of the electron fraction in the outflow. The mass flux carried by such outflows is found 
to be relatively large, on the order of 10~^Mq s~^. 

The main conclusion is that neutron rich outflows can be expelled from neutrino-cooled disks 
under certain conditions. However, those outflows are expected to be sub or at best mildly rela- 
tivistic, at least in cases where the central black hole is non rotating. Relativistic outflows can, in 
principle, be launched magnetically in the polar region, but those are, in general, proton rich. The 
above results may be altered if the central black hole is rapidly rotating. 

The slow, neutron rich winds discussed above may play an important role in the collimation 
of the central GRB-producing jet, and may also provide a contaminating baryon source (Levinson 
& Eichler 2003; McKinney 2006) for a central, baryon- free fireball that may be produced, e.g., 
via neutrino annihilation on magnetic field lines penetrating the horizon. Baryon pick-up by the 
central jet may have some interesting implications for fireball physics (e.g., Eichler &: Manis 2007). 

This work was supported by an ISF grant for the Israeli Center for High Energy Astrophysics 
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Fig. 1. — Flow parameters along a streamline for a pressure driven flow, assuming tanO — 0.2, surface 
magnetic field Bpo = 10^^ G, {B^/Bp)o = 0.1, and initial density po = 10^° gr cm~^. Each panel gives the 
slow magnetosonic Mach number Mgu (solid line), the dimensionless entropy per baryon s (dashed line), 
the normalized temperature and density (dotted lines), the electron fraction Y,. (dottcd-dashed line), and 
the free nucleon mass fraction X^ee, as functions of the normalized height above the disk z/Rq. 
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Fig. 3. — Same as fig. 1 for outflow from a turbulent disk. Initial values are indicated. The corresponding 
mass fluxes, as defined in eq. are M = 8 x lO^^, 2.7 x 10^°, 3.5 x lO^", 1.8 x 10^^ gr s~^ for the upper 

left, upper right, lower left and lower right panels, respectively. 



